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We present a measurement of the ratio of multijet cross sections in pp collisions at \fs = 1.96 TeV 
at the Fermilab Tevatron Collider. The measurement is based on a data set corresponding to an 
integrated luminosity of 0.7 fb -1 collected with the DO detector. The ratio of the inclusive three-jet 
to two-jet cross sections, R3/2, has been measured as a function of the jet transverse momenta. The 
data are compared to QCD predictions in different approximations. Popular tunes of the pythia 
event generator do not agree with the data, while SHERPA provides a reasonable description of the 
data. A perturbative QCD prediction in next-to-leading order in the strong coupling constant, 
corrected for non-perturbative effects, gives a good description of the data. 

PACS numbers: 13.87.Cs,12.88.Qk 



In hadron-hadron collisions, production rates of colli- 
mated sprays of hadrons, called jets, with large transverse 
momenta with respect to the beam axis (pt) are sensi- 
tive to both the dynamics of the fundamental interaction 
and to the partonic structure of the initial-state hadrons. 
The latter is usually parameterized in parton distribution 
functions (PDFs) of the hadrons. We observed large sen- 
sitivities to PDFs in our measurement of the differential 
three-jet cross section as a function of the three-jet in- 
variant mass [l|. Studies dedicated to the dynamics of 
the interaction are therefore preferably based on quan- 
tities which are minimally sensitive to the PDFs. Such 
quantities can be constructed as ratios of cross sections, 
for which the sensitivity to the PDFs is reduced. One 
class of such quantities is the ratio of multijet cross sec- 
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tions. The two-jet [2| , three-jet Q, and four-jet Q cross 
sections have been computed in perturbative Quantum 
Chromodynamics (pQCD) up to next-to-leading order 
(NLO) in the strong coupling constant a s . The ratio of 
the inclusive three-jet cross section to the inclusive two- 
jet cross section, -R3/2, provides a test of the correspond- 
ing NLO pQCD predictions. Previous measurements of 
R3/2 in processes with initial state hadrons, have been 
made in ep collisions at the HERA Collider at DESY [5| , 
and in hadron-hadron collisions at the SPS Collider at 
CERN 0, the Fermilab Tevatron Collider Q, and at the 
Large Hadron Collider at CERN [|, ||. 

This letter presents the first measurement of R3/2 in PP 
collisions at a center-of-mass energy of y/s = 1.96 TeV. 
The results are presented as a function of the highest 
jet pr in the event, prmax, for four minimum values of 
the second highest and (for three-jet events) third high- 
est jet pt, PTmin- The data sample, collected with the 
DO detector during 2004-2005 in Run II of the Fermilab 
Tevatron Collider, corresponds to an integrated luminos- 
ity of 0. 71b- 1 . 

Jets are defined by the Run II midpoint cone jet algo- 
rithm [l(| with a cone radius of i? C onc = 0.7. Rapidity is 
related to the polar scattering angle with respect to the 
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beam axis by y = 0.5 In [(1 + j3 cos0)/(l — (3 cos 8)], with 
j3 = \p\/E. The inclusive ra-jet event sample (for n = 2, 3) 
is defined by all events in which the n highest pt jets have 
Pt > PTmin and \y\ < 2.4. The separations in the plane of 
rapidity and azimuthal angle </>, Rjj — (Ay) 2 + (A0) 2 
between the n highest pt jets are required to be larger 
than twice the cone radius (Rjj > 2R conc ). The rapid- 
ity requirement restricts the jet phase space to the re- 
gion where jets are well-reconstructed in the DO detector 
and the energy calibration is known to 1.2-2.5% for jets 
with 50 < pr < 500 GeV. The separation requirement 
strongly reduces the phase space for which the n highest 
Pt jets had overlapping cones which were split during the 
overlap treatment of the jet algorithm. 

The ratio of inclusive three-jet to two-jet cross sections, 



-^3/2 (PTmax , PTmin) 



da 3 -}et(PTmin) / dp T 
d<T2-jct(pTmin) / dpr 



is less sensitive to experimental and theoretical un- 
certainties than the individual cross sections due 
to cancellations of correlated uncertainties. Here 
^3/2(Prmax, PTmin) is measured for p Tm \n requirements 
of 30, 50, 70, and 90 GeV in bins of pT max , in the in- 
terval 80 GeV < pTmax < 500 GeV, with the additional 
requirement of PTmax > PTmin + 30 GeV. This additional 
requirement ensures that there is sufficient phase space 
for the second and third jet so that corrections due to 
the experimental pt resolution remain small. Given the 
definitions above for the inclusive n-jet event samples, at 
each pTmax value the inclusive three-jet event sample is 
a subset of the inclusive two-jet event sample. Therefore 
^3/2(PTmaxj PTmin) represents the conditional probabil- 
ity for a two-jet event (at PTmax) to contain a third jet 

With p T > PTmin- 

A detailed description of the DO detector can be found 
in Ref. The event selection, jet reconstruction, jet 



energy and momentum correction in this measurement 
follow closely those used in our recent measurements of 
inclusive jet, two-jet, and three-jet distributions 0, [T2I — 
flH | . The primary tool for jet detection is the finely seg- 
mented uranium-liquid argon calorimeter that has almost 
complete solid angle coverage 1.7° < 9 < 178.3° [llj |. 
Events are triggered by a single high pt jet above a par- 
ticular threshold. In each PTmax bin, events are taken 
from a single trigger which is chosen such that the trig- 
ger efficiency is above 99% for two-jet and for three- 
jet events. Using triggers with different prescale val- 
ues results in integrated luminosities of 1.54 pb -1 for 
PTmax < 120 GeV, Upb- 1 for 120 < pr max < 140 GeV, 
73pb~ 1 for 140 < prmax < 175 GeV, 0.5 fb -1 for 175 < 
PTmax < 220 GeV, and 0.7 fib" 1 for p Tmax > 220 GeV. 

The position of the pp interaction, reconstructed us- 
ing a tracking system consisting of silicon microstrip de- 
tectors [13] and scintillating fiber tracker located inside 
a 2 T solenoidal magnet [ll| , is required to be within 
50 cm of the detector center along the beam direction. 



The jet four-momenta are corrected for the response of 
the calorimeter, the net energy flow through the jet cone, 
energy from event pile-up and multiple pp interactions, 
and for systematic shifts in rapidity due to detector ef- 



fects [15j. Cosmic ray backgrounds are suppressed by 



requirements on the missing transverse momentum in an 
event 15]. Requirements on characteristics of the shower 
shape are used to suppress the remaining background due 
to electrons, photons, and detector noise that mimic jets. 
The efficiency for these requirements is above 97.5%, and 
the fraction of background events is below 0.1% at all 

PTmax- 

The jet four-momenta reconstructed from calorimeter 
energy depositions are then corrected, on average, for the 
response of the calorimeter, the net energy flow through 
the jet cone, additional energy from previous beam cross- 
ings, and multiple pp interactions in the same event, but 

3 0. 



not for muons and neutrinos 



The absolute en- 
ergy calibration is determined from Z — > e + e~ events 
and the pt imbalance in 7 + jet events in the region 
\y\ < 0.4. The extension to larger rapidities is derived 
from dijet events using a similar data-driven method. In 
addition, corrections in the range 2-4% are applied that 
take into account the difference in calorimeter response 
due to the difference in the fractional contributions of 
quark and gluon-initiated jets in the dijet and the 7 + 
jet event samples. These corrections are determined us- 
ing jets simulated with the pythia event generator [l8[ 
that have been passed through a GEANT-based detector 
simulation The total corrections of the jet four- 

momenta vary between 50% and 20% for jet pt between 
50 and 400 GeV. An additional correction is ap plie d for 
systematic shifts in \y\ due to detector effects |l2l . Il5j |. 
These corrections adjust the reconstructed jet energy to 
the energy of the stable particles that enter the calorime- 
ter except for muons and neutrinos. 

The -R3/2 distributions are corrected for instrumental 
effects using a simulation of the DO detector response 
based on parameterizations of resolution effects in pt, 
the polar and azimuthal angles of jets, and jet recon- 
struction efficiencies. The parameterizations are deter- 
mined either from data or from a detailed simulation 
of the DO detector using GEANT. The parameterized 
simulation uses events generated with SHERPAvl.1.3 (2(| 
(including the tree-level matrix elements for two-jet, 
three-jet, and four-jet production) using default settings 
and MSTW2008LO PDFs and a sample of events, 
generated with PYTHlAv6.419 using tune QW [13] and 
CTEQ6.6 PDFs [13]. The events are subjected to the 
detector simulation and are reweighted such that their 
simulated distributions describe the differential two-jet 
and three-jet cross sections in the pt and rapidity of 
each of the three highest pt jets in the data. To min- 
imize migrations between PTmax bins due to resolution 
effects, we use the simulation to obtain a rescaling func- 
tion in PTmax that optimizes the correlation between the 
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reconstructed and true values. The rescaling function is 
applied to data and simulation. The bin sizes in pTmax 
are chosen to be much larger than the resolution. The 
bin purity after £>Tmax rescaling, defined as the fraction 
of all reconstructed events that are generated in the same 
bin, is above 50% for the two-jet and above 45% for the 
three-jet event samples. Bin efficiencies, defined as the 
fraction of all generated events that are reconstructed in 
the same bin, are above 55% for the two-jet and above 
45% for the three-jet event samples. 

We use the simulation to determine correction factors 
for the differential two-jet and three-jet cross sections in 
all PTmax bins, taking the average of SHERPA and PYTHIA. 
These include corrections for all instrumental effects, in- 
cluding the energies of unreconstructed muons and neu- 
trinos inside the jets. The total correction factors for the 
differential cross sections are between 0.92 and 1.0 for the 
two-jet and in the range 0.98-1.1 for the three-jet event 
samples. The correction factors for the ratio R3/2 are in 
the range 0.9-1.2. Over most of the range, the correc- 
tions from the two models agree within 3%. We take half 
the difference as an estimate of the model dependence 
of the correction, taking into account the correlations 
between the uncertainties for the two sets of correction 
factors. The corrected data are presented at the "particle 
level" (jets formed from stable particles after fragmenta- 
tion) as defined in Ref. [24j |. 

The measurement results are listed in Table U 

and Ref. [25J, and are displayed in Fig. Q] as a function of 
PTmax for different pTmin requirements. The ratio R3/2 
increases with increasing pTmax up to a maximum value 
and decreases for higher pr max values. The position and 
the height of the maximum depend on the pTmin require- 
ment (for the pTmin choices in this analysis, the maximum 
appears at PTmax values in the range 200-300 GeV). For 
a given value of PTmax, three-jet final states have on av- 
erage larger invariant masses than two-jet final states. 
Therefore the three-jet cross section approaches the kine- 
matic limit at lower prmax than the two-jet cross section, 
resulting in the decrease of R3/2 at large PTmax- The ini- 
tial increase of R3/2 with PTmax reflects the increasing 
phase space for three-jet final states, for a given PTmin 
requirement. For higher prmin requirements, the initial 
increase of R3/2 occurs at higher PTmax values, thereby 
shifting the position of the maximum. 

Theoretical calculations for R3/2 are computed as the 
product of NLO pQCD results and correction factors for 
non-perturbative effects. Predictions of NLO pQCD are 
obtained from NLOJET+- 1- j2(| using FASTNLO [27j] , Jets 
are reconstructed using the Fast Jet [28j implementation 
of the DO Run II midpoint cone jet algorithm. We use 
the two-loop approximation of the renormalization group 
equation for five quark flavors with a s (Mz) = 0.1180 



TABLE I: The ratio R3/2 measured as a function of PTmax 
for different pTmin requirements, along with statistical and 
systematic uncertainties. 



PTmax 

(GeV) 


PTmin 

(GeV) 


R3/2 




Stat, uncert. 
(percent) 


Syst. uncert 
(percent) 


80-100 


30 


1.816 


X 


10" 


1 


±0.7 


+5.6 


-5.5 


100-120 


30 


2.182 


X 


10" 


1 


±0.6 


±4.5 


-4.4 


120-140 


30 


2.370 


X 


10" 


1 


±0.5 


±3.7 


-3.7 


140-165 


30 


2.442 


X 


10" 


1 


±0.6 


±3.3 


-3.3 


165-190 


30 


2.464 


X 


IO- 


1 


±1.0 


±3.1 


-3.2 


190-220 


30 


2.421 


X 


D- 


1 


±0.6 


±3.1 


-3.1 


220-250 


30 


2.362 


X 


IC)" 


1 


±0.9 


±3.1 


-3.1 


250-285 


30 
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which is close to the world average value of 0.1184 29 1 

n, 



Results are computed using the MSTW2008NLO 
the CT10 HH, and the NNPDF2.1 [32] PDF sets. For 



consistency, we always use those PDFs which have been 
obtained for a s (Mz) — 0.1180. The renormalization and 
factorization scales [Ir and [If are set to hq = pTmax- 
The scale uncertainties are computed by varying and 
Hf independently between [io/2 and 2^to with the restric- 
tion that 0.5 < fiu/ Vf < 2.0. The uncertainties of the 
pQCD predictions due to the scale dependence are be- 
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FIG. 1: The measured -R3/2 results, compared to the predictions from NLO pQCD corrected for non-perturbative effects (top), 
and the ratio of data to theoretical predictions (bottom). The results are presented as a function of the highest jet pr, PTmax, 
for different pTmin requirements. The inner uncertainty bars represent the statistical uncertainties while the total uncertainty 
bars represent the quadratic sums of statistical and systematic uncertainties. (N.B. the inner uncertainty bars are within the 
markers for most data points.) 



tween -15% and +5%. 

The non-perturbative correction factors are the prod- 
ucts of hadronization and underlying event corrections. 
Both are estimated using pythia with tunes DW and 
AMBT1 [H]. Tune DW uses Q 2 -ordered parton showers 



TABLE II: The \ 2 values between theory and data for differ- 
ent pTmin requirements, for different choices of /j,r, hf- 



and an older underlying event model, while AMBT1 uses 
PT-ordered parton showers and a newer underlying event 
model. For each of these tunes, three event samples have 
been generated: parton shower level without an underly- 
ing event, particle level without an underlying event, and 
particle level with an underlying event. The hadroniza- 
tion corrections are estimated as the ratio of R3/2 at the 
particle level and at the parton level (from the partons 
at the end of the parton shower). Both are obtained 
without an underlying event. The underlying event cor- 
rection is the ratio of the particle level results with and 
without an underlying event. We use the average of the 
corrections obtained with tunes DW and AMBT1 as the 
central choice, and quote half the spread as the uncer- 
tainty. The total non-perturbative correction factors are 
in the range of 0.96-0.99 with uncertainties of less than 
1%. 

The theoretical predictions for the MSTW2008NLO 
PDF sets are overlaid on the data in Fig. [TJ The results 
for CT10 and NNPDFv2.1 PDFs (not shown) agree with 
those obtained for MSTW2008NLO to better than 0.1% 
for PTmax < 300 GeV, and are always better than 0.4%. 
Figure [TJ shows good agreement between the theoretical 
predictions and data. At the lowest PTmin value, the ra- 



PTmin 


Number of 


x 2 


for fj, R = pf 






data points 


PTmax /2 


PTmax 


2pTmax 


30 GeV 


12 


46.4 


21.7 


14.0 


50 GeV 


12 


12.4 


8.5 


9.1 


70 GeV 


11 


10.9 


9.6 


13.5 


90 GeV 


10 


13.3 


12.7 


14.4 



tio of data and theory decreases from +8% to —6% with 
increasing PTmax- For prmin = 50 and 70 GeV, the ratio 
of data to the theoretical predictions is consistent with 
unity over the entire range of prmax- For PTmin = 90 GeV 
the theoretical predictions are slightly higher than the 
experimental results, but still consistent within the the- 
oretical uncertainty. The agreement between theory and 
data is quantified by computing y 2 values for each choice 
of PTmin- The \ 2 definition takes into account all ex- 
perimental uncertainties and their correlations as well as 
uncertainties in the non-perturbative corrections and the 
PDFs. The y 2 values are listed in Table ITll for the differ- 
ent PTmin requirements and for three choices of /ir and 
Hf- For PTmin = 30 GeV, the x 2 value depends strongly 
on the scale, and agreement within the expectation of 
X 2 = V data ± v 72 • A^data (where V da ta is the number 
of data points) is obtained only for the largest scale of 
Hr,f = 2pTmax- For larger requirements of prmin = 50, 
70, and 90 GeV, the theoretical predictions agree with 
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FIG. 2: The measured R3/2 results, normalized to the predictions of the sherpa Monte Carlo event generator. The inner 
uncertainty bars represent the statistical uncertainties while the total uncertainty bars represent the quadratic sums of statistical 
and systematic uncertainties. Overlaid are the predictions from the pythia Monte Carlo event generator for four different tunes, 
also normalized to the sherpa predictions. 



the data for all three scales, and the central scale choice 
of fiR,F = PTmax always provides the lowest \ 2 - 

Predictions from different Monte Carlo event genera- 
tors are compared to the data in Fig. [2j In this Figure, 
the measured R3/2 results and the PYTHIA predictions 
for different tunes are divided by the predictions from 
sherpa which includes the tree-level matrix elements for 
two-, three-, and four-jet production, matched with a 
parton shower. The SHERPA predictions for R3/2 have a 
different prmax dependence and, independent of pTmin, 
they are approximately 20% lower (10% higher) than the 
data at low (high) pTmax- 

pythia includes only the two-jet matrix elements and 
a parton shower. The pythia predictions for the three- 
jet cross section therefore depend directly on the par- 
ton shower model and the corresponding parameter set- 
tings (tunes). The PYTHIA results have been obtained for- 
tunes DW, BW [H, A [13, AMBT1, S Global [H, Pe- 
rugia 2011, Perugia 2011 LO**, and Perugia 2011 Teva- 
tron [3^. The first three tunes use Q 2 -ordered parton 
showers and an older underlying event model, while the 
latter five use p^-ordered parton showers and a newer 
underlying event model. Tune DW was tuned to de- 
scribe the DO measurement of dijct azimuthal decorre- 
lations [37J, and tunes AMBT1, S Global and Perugia 
2011 were tuned to LHC data at -/s = 7TeV. The pre- 
dictions for tune Perugia 2011 agree within 1% with those 
for tunes Perugia 2011 LO** and Perugia 2011 Tevatron 
(the latter two are not shown in Fig. [2]). The predic- 
tions for tune A (not shown in Fig. [2]) are always above 
those for tune DW, and the predictions for tune S Global 
(not shown Fig. [5]) are 2-5% higher than those for tune 
AMBT1. Figure [5] shows that none of the studied pythia 



tunes describe the data; all predict a different pTmax de- 
pendence, and the discrepancies are strongly depending 
on pTmin- While a dedicated study of the sensitivity of 
the PYTHIA parameters is beyond the scope of this letter, 
these R3/2 data demonstrate the limitations of current 
pythia tunes and provide strong constraints for future 
parameter adjustments. 

In summary, we have presented the first measurement 
of the ratio R 3 / 2 of three-jet to two-jet cross sections in 
hadron-hadron collisions at a center of mass energy of 
y/s = 1.96 TeV. The ratio R3/2 is presented for PTmin re- 
quirements of 30, 50, 70, and 90 GeV, as a function of 
the highest jet pr, PTmax, in the range of 80-500 GeV. 
sherpa predicts a slightly different PTmax dependence, 
but it describes the data within approximately —10% to 
+20%. None of the pythia tunes DW, BW, A, AMBT1, 
S Global, and Perugia 2011 describe the data. The data 
are well described by the pQCD predictions at the next- 
to- leading order in the strong coupling constant a s , cor- 
rected for non-perturbative effects. 
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